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Abstract. We report a development of a multi-color simultaneous camera for the 188cm telescope at Okayama Astro- 
physical Observatory in Japan. The instrument, named MuSCAT, has a capability of 3-color simultaneous imaging in 
optical wavelength where CCDs are sensitive. MuSCAT is equipped with three 1024x 1024 pixel CCDs, which can be 
controlled independently. The three CCDs detect lights in P 2 (400-550 nm), (550-700 nm), and Zs ,2 (820-920 nm) 
bands using Astrodon Photometries Generation 2 Sloan filters. The field of view of MuSCAT is 6.1 x 6.1 aremin^ with 
the pixel scale of 0.358 arcsec per pixel. The principal purpose of MuSCAT is to perform high precision multi-color 
transit photometry. For the purpose, MuSCAT has a capability of self autoguiding which enables to fix positions of 
stellar images within ^1 pix. We demonstrate relative photometric precisions of 0.101%, 0.074%, and 0.076% in g 2 , 
?’ 2 , and Zs 2 bands, respectively, for GJ436 (magnitudes in p'=l 1.81, r'=10.08, and z'=8.66) with 30 s exposures. The 
achieved precisions meet our objective, and the instrument is ready for operation. 
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1 Introduction 

Transiting planets, which transit in front of their host stars, are especially important research ob¬ 
jects among exoplanets, as when combined with radial velocity measurements, they can provide us 
various information about the nature of exoplanets such as the mass, radius, density, orbital obliq¬ 
uity, and atmosphere. Most of transiting planets have been discovered by transit surveys which 
monitor the brightness of hundreds of thousands of stars. Several groups have worked or been 
working on ground-based transit surveys,and CoRoT^ and Kepler^ have performed space-based 
transit surveys. Moreover, the second epoch mission of Kepler, namely K2,® is now ongoing, and 
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next generation spaee missions TESS’ and PLATO^ have been approved to be launehed around 
2017 and 2024, respeetively. 

The biggest problem for transit surveys, espeeially for ground-based ones, is eommingling 
of false positives with eandidates of transiting planets. The eause of false positives of transiting 
planets are eelipsing binaries, as they mimic planetary transits by grazing or by being mixed in 
another bright star’s brightness. Thus it is important to conduct follow-up observations after transit 
surveys in order to validate candidates as true planets by eliminating false positives. 

High precision multi-color transit photometry is known to be useful for such follow-up obser¬ 
vations to discriminate whether transit-like dimming is caused by a true planet or by an eclipsing 
binary.^ This is because a true planet is almost dark in all wavelength, while an intervening body of 
an eclipsing binary is bright itself and its brightness changes significantly with wavelength. Thus 
false positives caused by eclipsing binaries can be spotted by observing significant wavelength 
dependence in transit depths. On the other hand, transit depths of a true planet also have wave¬ 
length dependence. Most of the wavelength dependence comes from the stellar limb-darkening, 
but the apparent planetary radius also has weak wavelength dependence which reflects the nature 
of its planetary atmosphere. High precision multi-color transit photometry is known to be useful 
to measure the weak wavelength dependence in transit depths to study atmospheres of transiting 
planets. This kind of study is known as transmission spectroscopy, and numbers of multi-color 
transit observations for this purpose have been reported so far.^° 

Multi-color simultaneous cameras are very fruitful for the studies described above for two 
reasons. First, multi-color simultaneous cameras provide not only higher efficiency but also more 
feasibility to accomplish aimed studies than single-color cameras, since observable transits from 
a specific ground site are very limited. Second, simultaneity of multi-color transit photometry is 
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important to avoid systematic differences of transit depths due to luminosity change in host stars 
possibly caused by existence of starspots, plages, stellar activity, and so on. For the reasons, multi¬ 
color simultaneous cameras such as GROND,^°’'^ BUSCA,^^ ULTRACAM,^^ SIRIUS,'"^’and 
MITSuME^^’^^ have been actively used for transit observations. 

Considering the fact that more interesting transiting planets will be discovered in the near fu¬ 
ture by advanced ground-based surveys, and also space-based surveys, like K2, TESS, and PEATO, 
developments of new multi-color simultaneous cameras are highly desired. We here report a devel¬ 
opment of such an astronomical instrument named MuSCAT (Multi-color Simultaneous Camera 
for studying Atmospheres of Transiting planets), which is now installed on the 188cm telescope at 
Okayama Astrophysical Observatory (OAO) in Japan. 

The rest of this paper is organized as follows. We first describe designs of the optical system 
of Muscat and its components (Sec. [2l), and introduce the control system of MuSCAT (Sec. [3]). 
We then report characteristics and performances of MuSCAT shown in engineering observations 
(Sec. IH). We discuss on some capabilities for future upgrade of MuSCAT (Sec. [5]), and finally 
summarize this paper (Sec. [b]). 

2 Optical System 

2 .1 Scientific Requirements and Design Policies 

We have designed MuSCAT considering the following conditions. As we plan to use MuSCAT 
for validations of transiting planets discovered by transit surveys, at least 2 colors are necessary to 
discriminate eclipsing binaries from transiting planets. Considering the cost and available research 
grants, we adopt a design for a 3-color simultaneous camera with the 3 colors in optical wavelength 
where CCDs are sensitive. Eor transit observations, it is important to obtain good comparison stars 
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Fig 1 Layouts of the optical system of MuSCAT; (a) in 3D and (b) in 2D. Meanings of the labels in the panels are as 
follows. A : Telescope focal plane (126mml26mm). B : F conversion optics, constructed in order with S-BSL7 lens, 
45deg mirror, CaF2 lens, S-LAL8 lens, CaF2 lens, and finally BK7 lens after a bandpass filter in front of a detector. C 
: Offner relay. D : Dichroic mirrors. E : Bandpass filters. F : Detectors. The vertical dashed lines in panel (b) indicate 
the same positions, but the bottom shows a view from the bottom. 


in the field of view (FOV) to achieve high precision transit photometry. For the reason, we have 
designed the MuSCAT FOV as wide as possible for the 188cm telescope. We have also took care 
of the throughput (TP) of the instrument to achieve high photometric precision. To achieve higher 
sensitivity to the utmost extent, we have carefully selected and designed the MuSCAT optical 
system including astronomical bandpass filters, dichroic mirrors, and CCDs. 


2.2 Optical Design 

Layouts of the optical system of MuSCAT are shown in Fig. [T] MuSCAT adopts a 45° plane mirror 
and an Offner relay, which consist of Si02-protected aluminum mirrors, and inserts 2 dichroic 
mirrors in the light path to take images of 3 bands in optical wavelength simultaneously. MuSCAT 
maintains space for another dichroic mirror around the Cassegrain focus so as to accommodate 
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Fig 2 A picture of MuSCAT installed on the 188cm telescope. 


near-infrared (NIR) channels as a future upgrade capability. For high precision transit photometry, 
it is desired to have the widest possible FOV to get good comparison stars. The Offner relay enables 
us to distribute F conversion lenses around the relay. The lenses contribute to achieve wider FOV 
by converting F-number from FI8 at the telescope to F5.5 before the relay and further converting 
to F4.0 just before CCD cameras after the relay. Aberration correction is performed by 4 lenses 
before and 1 lens after the Offner relay, and finally accomplished by a plane-convex lens just before 
CCD cameras. Anti-reflection coating for optical wavelength is applied to the lenses. The lenses 
are designed to correct aberration of the whole system including both the 188cm telescope and 
Muscat. Thereby the optical system of MuSCAT provides a good imaging quality throughout the 
FOV. Wavelength divisions are performed by 2 dichroic mirrors after the Offner relay. The dichroic 
mirrors are wedge-shaped to reduce astigmatism. Astronomical bandpass filters are inserted just 
before the last lenses. Lights of astronomical objects are divided into 3 colors and detected by 3 
CCD cameras. A picture of the actual MuSCAT installed on the 188cm telescope is presented in 
Fig.E 
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2.3 Dichroic Mirrors and Bandpass Filters 


The 2 dichroic mirrors are manufactured by Asahi Spectra Co.,Ltd. The size of the first dichroic 
mirror (DM1) is 113 mm by 108 mm, and the depth of 11.8 mm with the wedge angle of 7 min 
52 sec. The size of the second dichroic mirror (DM2) is 90 mm by 88 mm, and the depth of 
9.9 mm with the wedge angle of 12 min 33 sec. Anti-reflection coating is processed on the back 
sides of the dichroic mirrors so that the DMs transmit remaining lights (namely, not reflected ones) 
almost completely. Fig. [3]plots wavelength dependence of reflectance of DM1 and DM2 measured 
by a gonio-spectrophotometer. The reflectance of both DMs is almost flat across the relevant 
wavelength for the 3 bands. We note that there are wiggles in reflectance around 350 nm for both 
DMs and around 525 nm for the DM2, but those wiggles do not affect to the current 3 bands. 

The DM1 is inserted with an incident angle of 16 deg and it reflects most of light in 400-550 
nm and transmits remaining light longward of 550 nm. Reflected light from the DM1 is detected 
by the eh 1 CCD camera. Then the DM2 is inserted with an incident angle of 22.5 deg and it 
reflects most of light in 550-700 nm and transmits remaining light longward of 700 nm. Reflected 
light from the DM2 is detected by the eh 2 CCD camera, and transmitted light from the DM2 is 
detected by the eh 3 CCD camera. 

We adopt § 2 , r^, Zs ^2 band filters of Astrodon Photometries Generation 2 Sloan filters. The size 
of those filters are 50 mm by 50 mm. Fig. 0] plots the wavelength dependence of transmittance of 
the bandpass filters, measured by a spectrophotometer SHIMADZU UV-3100PC. The reflectance 
and transmittance of DM1 and DM2 are optimized for those 3 bands. 

We note that changes in the reflectance and transmittance of DMs and bandpass filters are 
negligible in normal operations unless dew condensations occur. 
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Fig 3 Wavelength dependence of reflectance of DM1 and DM2. 

2.4 Cameras and Detectors 

Muscat equips 3 CCD cameras manufactured by Princeton Instruments. The first one is a 
PIXIS: 1024B model, used as the eh 2 CCD camera in band. The other two are PIXIS: 1024B_eXcelon 
model cameras, used as the eh 1 {g 2 band) and eh 3 {Zs ,2 band) CCD cameras. Each CCD camera 
equips back-illuminated grade 1 CCD chip (e2v CCD47-10 for PIXIS: 1024 and Princeton Instru¬ 
ments’ proprietary CCD for PIXIS: 1024B_eXcelon) with Ikx Ik (1024x 1024) pixels. It is noted 
that the CCDs of PIXIS: 1024B_eXcelon are specially-processed to suppress the etaloning (fring¬ 
ing) that occurs in standard back-illuminated CCDs. Quantum efficiencies (QE) of PIXIS: 1024B 
and PIXIS: 1024B_eXcelon are plotted in Pig. [51 Special anti-reflection coating (BBAR coating) 
is applied to the vacuum windows of CCD cameras and its transmittance is shown in Pig. The 
data of the QEs and the transmittance of the BBAR coating are provided by the manufacturer. We 
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Fig 4 Wavelength dependence of transmittance of bandpass filters, 
summarize nominal specifications of the CCD cameras in Table [B We also present actual mea¬ 
sured values of gains, full well, and read noise as well as an upper limit of dark current based on 
data taken during engineering observations (Sec. S]) in Table [B 

2.5 Total Throughput 

Efficiencies of the F conversion optics and the Offner relay are roughly estimated as 60% in 
band, 61% in band, and 50% in Zgg band. Based on the transmittance and reflectance of DMs, 
filters, BEAR coating, and QE of CCDs in the previous subsections, we calculate expected total 
throughput of MuSCAT. We plot the wavelength dependence of the expected total throughput of 
Muscat in Fig. Hand present machine-readable values in Table[2l A comparison of the expected 
total throughput with a measured one is presented in section 4.3. 
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Fig 5 Wavelength dependence of quantum efficiency of CCDs. 

2.6 Field of View 

Muscat has 6.1 x6.1 arcmin^ FOV with the Ikx Ik CCDs introduced above. The pixel scale for 
each CCD is ~0.358 arcsec per pixel (see Fig[T8]for details). Centers of FOV of 3 CCD cameras 
are adjusted within 10 pixels in both x and y directions, and relative rotations of the position angle 
fit within 1 deg. 

As a future upgrade capability, MuSCAT can replace CCD cameras with 2kx2k (2048x2048 
pixels) CCD cameras PIXIS: 2048B and PIXIS: 2048B_eXcelon. With such an upgrade, the FOV 
of Muscat will get larger to 12.7 x 12.7 arcmin^. In that case small (0-4%) vignetting depending 
on wavelength would occur in the region outside the diameter of 14.4 arcmin from center of FOV, 
while there is no vignetting for Ikx Ik CCDs. Fig. [8] indicates the explained potential FOV of 
Muscat. 
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Fig 6 Wavelength dependence of transmittance of BEAR coating. 

2.7 Spot Diagram 

Simulated imaging performances of the MuSCAT optical system are shown in spot diagrams in 
Figs. |9] and [TOl Fig. |9] plots spot diagrams for on-focus cases which indicate spot radius of 
all wavelength are well less than 1 arcsec throughout the current FOV and the potential FOV. 
This imaging performance is thus sufficient for the Okayama Astrophysical Observatory where 
the typical seeing is about 1.5 arcsec. Fig. [TOl shows spot diagrams for defocused cases where the 
secondary mirror is shifted by 1.5 mm, which makes spot radius expand to about 4 arcsec. It is well 
known that defocusing is very useful for high precision transit photometry for isolated sources,*^ 
and thus those cases are more realistic for transit observations. The panels imply that images are 
almost circular throughout the FOV and suitable for aperture photometry. 
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Fig 7 Wavelength dependence of expected total efficiencies of MuSCAT in g 2 , and 2^,2 bands. 



Fig 8 A schematic view of the potential held of view of MuSCAT. 
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(a) 


(b) 


(c) 


Fig 9 Spot diagram of the CCD cameras for on-focus cases for (a) ch 1, (b) ch 2, and (c) ch 3. The 5x5 cells represent 
the FOV of 2kx2k CCD case, while the inner 3x3 cells do a case for Ikx Ik CCD. The size of each cell corresponds 
to 2x2 arcsec. Colors indicate simulated images of the shortest (blue), mid (green), and longest (red) wavelength in 
each channel. Specifically, 400, 470, 550 nm for ch 1, 550, 630, 700 nm for ch 2, and 700, 800, 950 nm for ch 3. 
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Fig 10 Same as Fig. |9l but for defocused cases. The size of each cell corresponds to 10x10 arcsec. 
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3 Control System 


3.1 System Structure 

Fig. [TT] shows a flow chart of the MuSCAT control system. The system consists of 2 personal 
computers: one (PCI) is directly mounted on the instrument and the other (PC2) is located in an 
observation room in the observatory. We use Linux/CentOS 5 (32bit) for PCI and Linux/CentOS 7 
(64bit) for PC2. For PCI, a 32bit operating system is required to customize observation commands 
using 32bit library for the CCD cameras. Those operating systems may be subject to change in the 
future. 

Each of the CCD cameras is connected to PCI with a USB cable. Each camera is independently 
controlled by PCI using an interface software which is developed based on PVCAM and CEITSIO 
C libraries. The images taken by these cameras are saved in a 16-bit PITS format along with PITS 
header information including the time information, telescope status, dome status, and weather 
information. Por the time information, PCI refers to the internal clock which is synchronized to 
Network Time Protocol servers via the internet. The typical time offset between the local machine 
and the NTP servers is less than 1 milliseconds. Thus total systematic uncertainties on the times 
of exposures, including the time lag between the internal time acquisition and actual exposure, 
are well within 1 second, which is negligible for even time-critical sciences such as transit timing 
measurements. Por the other information than the time, PCI refers to a telescope-control PC, 
which gathers all the up-to-the-second information. 

PC2 provides a user interface. Observers can send an observing command to each camera 
from this PC specifying the exposure time, number of images, gain and readout speed settings, 
and so on. We note that one can select optimal settings for each camera. Namely, one can set a 
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• Data storage 


Fig 11 A flow chart of MuSCAT control system. 

different exposure time and a different readout speed to each camera. This is an important feature 
of Muscat that we can use an optimal exposure time for each band independently. The observed 
images are instantaneously sent to PC2 to be displayed in SAOImage DS9. PC2 also has a function 
of self autoguiding (see Section [T2l) . and has a data storage. 


3.2 Self Autoguiding 

To achieve a 0.1% level high precision relative photometry, it is essential to receive stellar fluxes 
by the same pixels during a set of observations to mitigate the incompleteness of the pixel-to-pixel 
sensitivity correction. Because the tracking of the 188cm telescope is not perfect, an autoguiding 
system is critical to keep the stars on the same positions on the detectors. However, MuSCAT has 
no guide camera to capture a guide star located in the surrounding area of the FOV of the science 
cameras. We therefore have developed a self autoguiding system which uses scientific images for 
guiding.Specifically, the stellar centroids of several bright stars are measured on one of the 
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three band images soon after the image is obtained. Then the mean stellar displacement on the 
latest image relative to a reference image is calculated to feed back to the telescope. All the above 
processes are done within a few seconds after the end of exposure. We note that the CCD channel 
to be used for the stellar centroid calculation is selectable. We usually set the exposure time for the 
guiding channel to 30-60 seconds such that the feedback result will be well reflected in the next 
image and that the guiding frequency is high enough. As shown in Section I4.4[ the autoguiding 
system can stabili z e the stellar centroid positions within ~1 pixel for bright stars (magnitude less 
than about 12). This is important to reduce systematic errors caused by incompleteness of flat- 
fielding. We find no large difference in the autoguiding performance between in-focus and out-of- 
focus observations, meaning that the guiding performance is not limited by the degree of defocus 
but limited by the change of stellar PSF shapes due to seeing variation and the mechanical accuracy 
of telescope driving (tracking and fine moving). 

4 Results of Engineering Observations 

We conducted the first light engineering observation on the night of December 24, 2014, and 
further engineering observations on the nights of March 2-4 and April 3-5, 2015. We examine the 
performance of MuSCAT using data taken on the nights. We summarize results of the engineering 
observations below. 

4.1 Detector Characteristics including Bias, Flat, and Linearity 

We have took hundreds of bias, flat, and linearity test frames during engineering observations in 
order to learn detector characteristics of MuSCAT. We have observed lights of a filament lamp 
projected onto a matte whiteboard on a wall of the 188cm telescope dome for flat frames and 
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Fig 12 Median flat images of (a) band (with low etaloning process) and (b) band (standard CCD). 


linearity test frames. We do not find any strange features or signifieant bad pixel regions in bias 
and flat frames, exeept a unique feature of low etaloning proeess in and Zs ^2 band flat images. 
For referenee, we present median flat images in 513 (median of 433 flat frames) and (median of 
490 flat frames) bands in Fig. [121 We note that the unique feature has a very good repeatability 
with the fractional fluctuation of much less than 0.1%, which has little impact on the photometric 
precision required for our purposes. We have also exposed a very bright star of 1/ = 6 to check for 
the existence of image persistence on the CCDs and ghost patterns due to reflection by the lenses. 
We found no apparent features that can affect the photometry. We have derived full well values, 
gains, readout noises for the data. The values are presented in Table HI 

We have also tested the linearity of MuSCAT CCDs for each readout speed and each gain. Our 
method is based on a previous study for the CCDs of High Dispersion Spectrograph (HDS) of the 
Subaru telescope.*^ First, we have created linearity test frames which have gradational counts on 
the CCDs, by opening only a half of the tertiary mirror cover and inserting a black plate into the 
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Fig 13 An example of a linearity test frame in p 2 band. Counts of the dark region are nearly bias level and those in 
the brightest region are saturated. 
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Fig 14 An example of a linearity test result. The horizontal axis (f(8s) [e-]) indicates electron counts in 8 seconds. 
The vertical axis shows the flux ratio parameter which is defined as f(8s) - 2 * f(4s) [e-]. 
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light path in front of MuSCAT. Fig. [13] shows an example of a linearity test frame. Second, we 
monitor counts of the filament lamp until the filament lamp is stabilized. We note that it takes about 
2 hours until counts become nearly-unchanged. We then start linearity test exposures as follows. 
We first determine an exposure time for each CCD which gives counts from bias level to saturation 
level gradationally on the CCDs. We define frames with the above exposure time as “A” frames 
and ones with a half of the exposure time as “B” frames. We then take A and B frames alternately 
until obtaining 20 frames each. We have repeated such exposures for each gain and each readout 
speed, namely for the gain modes of 1, 2, 4 e“/ADU, and for the readout speed of lOOkHz and 
2MHz. Subsequently we subtract a median bias frame for each gain and each readout speed. We 
then make a new frame which computes photon counts of each pixel in an A frame minus twice 
of photon counts for the same pixel in a B frame using adjacent A and B frames (39 pairs in total 
for each gain and each readout speed). We define those frames as “C” frames (namely, C = A — 
2 X B for each pixel). To visualize the linearity of the CCDs, we plot electron counts (namely, 
photon counts x gain) of pixels in A frames as X-axis and electron counts of the same pixels in 
C frames as Y-axis. An example of such a figure is shown in Fig. [T4| We finally fit the plotted 
data with a linear function (Y = aX) using the data up to X=64000, and the best fit linearity slopes 
are summarized in Table (3] Based on the above test, we have confirmed that MuSCAT CCDs have 
a good linearity within ~0.21% at a maximum up to the saturation level. The result means that 
the effect of non-linearity is well negligible even for high precision transit photometry if counts of 
stars do not change drastically during observations. In the case we need to correct non-linearity, 
we will use those data for non-linearity corrections. 
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Fig 15 PSF on the detectors extracted from M67 images in (a) 53 band, (b) band, and (c) 2^,2 band. 
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Fig 16 Same as Fig. [15] but ellipticity maps on the detectors. 
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Fig 17 Same as Fig. [15] but FWHM maps on the detectors. 


19 











































+04’ 

+0.V 
+ 02 ' 

+or 
+00' 

■or 
-02’ 

-03' 

-04' 

Fig 18 Same as Fig. [15] but pixel scale distortion maps on the detectors. 

4.2 PSF and Distortion 

We obtained images of the open cluster M67 with MuSCAT during the first light observation. 
On this night, the sky condition was photometric and there was no moon. The airmass of M67 
was ~1.2. Using the M67 images, we have derived point spread function (PSF) and distortion 
of Muscat images using IRAF, PSFEx,^° SExtractor,^^ and SCAMP.^^ The extracted PSF, its 
ellipticity and FWHM of MuSCAT CCDs are shown in Figs. [151 [Ml and [T^ We note that we 
have used ”PIXEL_AUTO” for the parameter of the PSE model in PSEEx. The software estimates 
representative PSE for each grid using stars inside the grid. We also note that the seeing in g 2 , r^, 
and Zs ^2 was 2.1", 1.8", and 1.7", respectively, which was slightly worse than the typical seeing at 
the site of ~1.5" in optical bands. Thus imaging quality was not limited by MuSCAT itself but 
by the seeing. We have confirmed that the PSE is nearly circular throughout the EOV and that 
MuSCAT does not have unexpected large aberration or imaging problems. 

We have also derived distortion maps of images as differences in the pixel scale on the CCDs 
using one of output options of SCAMP.^^ The derived pixel scale distortion maps are presented in 
Eig. mi The figure indicates that the pixel scale distortion is limited within about 0.3%, which is 
negligible for standard aperture photometry. 
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4.3 Sensitivity and Efficiency 


We also estimate limiting magnitudes for the g 2 , r 2 , and Zs ,2 bands using the images of M67. The 
measured sky brightness in g!^, and Zs ,2 were 19.9 mag aresee”^, 19.5 mag aresee”^, and 18.7 
mag arcsee”^, respeetively. For each band, 10 x 60 s images were obtained with dithering. We 
note that we used high speed readout mode (2 MHz). We apply bias-fiat correction and stellar 
position alignment to the data and stacked them into a single image for each band. We conduct 
photometry for ~100 stars on each stacked image by using the DoPHOT package,^"* which per¬ 
forms an analytical PSF fitting. The measured instrumental magnitudes of these stars were then 
compared to the SDSS 9 catalog^^ for photometric calibration. We note that we here ignore color 
terms and simply approximate that the ^2’ ^s ,2 bands are identical to the g', r', and z' bands, 

respectively. Limiting magnitudes with 10-min exposure are estimated as the signal-to-noise (S/N) 
ratio reaches 10, yielding g[^.^ = 21.7, r[:^^ = 21.7, and ziim = 19.8. We note that zu^ is affected 
by higher readout noises in Zs ^2 band (see Table [1]) as we used the high speed readout mode. For 
the S/N calculation, we simply adopt the photometric errors returned by DoPHOT. We show a plot 
for SDSS magnitudes v.s. photometric errors in Fig. [T^ 

In addition, total throughput including the airmass, 188cm telescope, and MuSCAT is measured 
by the same data. We estimate total throughput as follows. First, we measure the zero-point mag¬ 
nitudes on the stacked M67 images as ZP(5'0=28.63, ZP(r2)=28.71, and ZP(2;s,2)=27.08, which 
correspond to 10 electrons for all bands. Next, we estimate the expected incident flux coming from 
an astronomical object with the above magnitudes into the effective area of the primary mirror of 
the 188cm telescope. Finally, comparing the expected flux with the detected one (10 electrons), 
we estimate the total throughput in 5(3, r^, and z* 2 bands as 20%, 28%, and 13%, respectively. 
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Mag vs Error for M67, Exp=60s x 10 



Fig 19 Relations of SDSS magnitudes and photometric errors for M67 stars taken with MuSCAT. A dashed line in 
magenta, a dotted line in orange, a dashed-dotted line in cyan, and a solid line in gray represent stellar photon noise, 
readout noise, sky photon noise, and total (theoretical) noise, respectively. The solid black line is the best-fit fifth-order 
function, just for visibility. 

The actual measured values are almost the same with expected values from the airmass, 188cm 
telescope, and MuSCAT, as shown in Table IH 

4.4 High Precision Photometry 

In order to check for the photometric performance of MuSCAT, we observed two stars; GJ 436^^’^^ 
and WASP-12,^^ both hosting a transiting planet. Observations were carried out on 2015 March 
2 UT, during out-of-transit phases for both targets. There was no cloud, but the sky level was 
relatively high due to a bright waxing moon with the age (lunar phase) of 11.5. We used the low 
speed readout mode. GJ 436 is a nearby M3.5 dwarf with 5^'=11.81, r'=10.08, and ^'=8.66 in 
the SDSS magnitudes. GJ 436 was observed during 14:50-15:57 UT (airmass=1.05-1.01) with 
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the exposure time of 30 s for all bands. WASP-12 is a GO dwarf with i?=12.1, R=\\.2 (from 
the NOMAD 1.0 eatalog^^), and z'=llAl?^ WASP-12 was observed during 13:37-14:39 UT 
(airmass=l .21-1.46) with the exposure time of 30 s (g'2 and r^) or 60 s (^s,2)- The exposure time and 
the number of exposures for each target, each filter is summarized in Table[51 For each observation, 
the FOV was adjusted so that several similar-brightness stars were simultaneously imaged. The 
stellar image was defocused such that the FWHM of PSF became ~33-38 pixels and ~31-36 
pixels for GJ436 and WASP-12, respectively. The self autoguiding system (see Sec. 13.21) was 
activated by using 5^2-band images for both targets. The stellar centroid changes in band during 
the observations of GJ 436 are displayed in Fig. [20l showing that the stellar positions were quite 
stable with the dispersion not exceeding ~1 pixel. 

The observed data were reduced by using a customized aperture-photometry pipeline.^' The 
applied aperture radius, the number of comparison stars used for the relative photometry, and the 
unnormalized flux ratio of the target star and the ensemble of the comparison stars are summarized 
in Table [51 Note that the applied aperture radius was determined such that the dispersion of the 
resultant light curve was minimum. We show the resultant normalized light curves of GJ436 and 
WASP-12 in Figs. |2T] and |22l respectively. The black dashed line indicates the best-fit linear 
function. Photometric precisions, which we define as the root-mean-square (rms) of the residual 
light curve from the liner fit, achieve 0.101%, 0.074%, and 0.076% in the g 2 , r^, and 2 bands for 
GJ 436, while those for WASP-12 are 0.16%, 0.16%, and 0.15%, respectively. These rms values 
are listed in Table [51 

To see how the photometric performance of MuSCAT has been achieved, we calculate the error 
budgets for these observations as shown in Tabled In the table, ^target, c’^comp, and CTgky indicate 
photon noises arising from the target-star flux, comparison-star flux, and sky-background flux. 
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respectively, calculated assuming the Poisson (photon) noise. (Tread is the noise contributed from 
the read-out noise listed in Table [IJ crgcin is the scintillation noise, for which we apply the following 
equation. 


Tscin = 0.064D-2/3(sec (1) 

where D is the diameter of the primary mirror of the telescope in cm, Z is the zenith distance, 
h = 372 m is the height above sea level of the observatory, ho = 8000 m, and T is the exposure time 
in seconds.All the remaining (unknown, or difficult to assess) components of the photometric 
error are treated as (Tunknown, which is calculated as 


^unknown 


rms^ 


a. 


target 


O'. 


comp 


a. 


sky 


cr: 


read 


a. 


scin’ 


( 2 ) 


where rms is the same as that listed in Table[5l Among these noise sources, all but (Tunknown are ba¬ 
sically unavoidable. Possible causes of (Tunknown can be the difference of atmospheric transparency 
between toward the target star and toward the comparison stars, the modulation of scintillation 
noise, the incompleteness of flat-field correction, and so on. We indeed find that (Tunknown is the 
major noise source in some of the light curves, but is still limited in degree of about 30-40% in 
rms^, meaning that (Tunknown is not a very limiting factor for the photometric precision. In other 
words, the most part of the photometric precision (> 60% in rms^ for all three bands) can be 
explained by the theoretical noise models. We therefore consider that the expected photometric 
performance of MuSCAT has been well achieved. 

We also note about the time-correlated noise (so-called the ’’red” noise) in the observed data. 
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r' 2 —band 



Fig 20 Time series centroid positions of bright stars in the field of view for GJ 436. 

For high-precision photometry such as transit observations, a treatment of the red noise would be 
very important.^"* We calculate a red-noise factor, which is the ratio of measured rms in binned data 
to the one expected from the rms in unbinned data, for our observations. We find 1.3 in average, 
which is a typical value for ground-based transit observations. Although one hour observations 
are not sufficient to evaluate the red noise in detail, we consider that the level of the red noise of 
Muscat is similar to other ground-based instruments. We will thus take into account the red noise 
for future science observations with MuSCAT. 

5 Upgrading and Transferring Capability 

Although the current MuSCAT is ready for operations, it still has upgrading capabilities. First, 
Muscat can be upgraded in terms of FOV by replacing the three Ikxlk CCD cameras, which 
give 6.1x6.1 arcmin^ FOV, with 2kx2k CCD cameras, which will provide 12.7x12.7 arcmin^ 
FOV. Such a wider FOV would be desirable to find good comparison stars especially for very 
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GJ436 



Fig 21 Light curves of GJ 436. 

bright targets which will be discovered by the TESS mission. Second, MuSCAT has space for 
another dichroic mirror to add NIR channels (see Sec. 12.21) . Additional NIR channels enable us 
to take images from optical to NIR simultaneously like GROND. Such a capability will enhance 
scientific merits for transmission spectroscopy in the light of efficiency and simultaneity. 

Muscat also has a transferring capability. The current instrument is optimized for the 188cm 
telescope at Okayama Astrophysical Observatory whose F-number is FI8, but MuSCAT can be 
transferred to or can make a copy of itself for other telescopes by replacing F conversion lenses. 


6 Summary 

We have developed a new astronomical instrument MuSCAT for the 188cm telescope at Okayama 
Astrophysical Observatory in Japan. MuSCAT has a capability of 3-color simultaneous imaging in 
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WASP-12 



Fig 22 Same as Fig. 121] but for WASP-12. 

5^2 (400-550 nm), r'^ (550-700 nm), and 2 (820-920 nm) bands with three Ikx Ik pixel CCDs. 
The field of view of MuSCAT is 6.1 x 6.1 aremin^ with the pixel seale of 0.358 aresee per pixel. 

One of the prime aims of MuSCAT is to eonfirm whether eandidates of transiting planets dis¬ 
covered by transit surveys, including such as K2, TESS, and PLATO, are truly planets or false 
positives due to eclipsing binaries. Another prime aim of MuSCAT is to measure the wavelength 
dependence of transit depths in visible bands, providing rough information about exoplanetary at¬ 
mospheres: such as the feature of the Rayleigh scattering by hydrogen dominated atmospheres, the 
feature of the Mie scattering by hazy atmospheres, or the flat feature of cloudy atmospheres. 

The capability of multi-color simultaneous transit photometry is well suitable for those aims. 
Since MuSCAT can achieve 0.1% photometric precision with 30 s exposure for stars brighter than 
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~10 mag as shown in Sec. I4.4[ MuSCAT will work effectively for the purposes especially for 
bright TESS transiting planets. In addition, MuSCAT would be also useful for follow-up ob¬ 
servations of supemovae and gamma-ray bursts, and monitoring variable stars, and so on. The 
instrument is ready for operation at Okayama Astrophysical Observatory. 
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Table 1 Summary of basic information of cameras and CCDs 


Camera name 

PIXIS 1024B / PIXIS 1024B eXeelon 

CCD format 

Pixel size 

Imaging area 

Data interfaee 

Dimensions 

Weight 

1024 X 1024 pixels 

13 X 13 /rm 

13.3 X 13.3 mm 

USB2.0 

16.59 em x 11.81 em x 11.38 em (L x W x H) 
2.27 kg 

Seleetable gain modes 
mode 1 

Measured gains mode 2 

mode 4 

Dark eurrent @ —70°C 

l,2,4e-/ADU 

1.00 (eh 1), 1.02 (eh 2), 1.00 (eh 3) e’/ADU 

1.98 (eh 1), 2.02 (eh 2), 2.00 (eh 3) e'/ADU 

3.96 (eh 1), 4.05 (eh 2), 3.97 (eh 3) e'/ADU 
<0.0015 e“s“^pix“^ (measured) 

0.0004 e“s“^pix“^ (eatalog typieal value) 

Full well: Single pixel 

Seleetable readout speed (readout time) 
Measured read noise: @ 100 kHz 

@2 MHz 

~70000 e- 

100 kHz (10.0 s), 2 MHz (0.58 s) 

3.9 (eh 1), 3.8 (eh 2), 4.2 (eh 3) e" 

11 (eh 1), 12 (eh 2), 24-27'^ (eh 3) e" 


An additional pattern (not random) noise is seen only in the 2MHz mode of the ch 3 CCD. This is scheduled to be repaired (the repair has been 

completed at the time of publication). 
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Table 2 Wavelength dependence of expected total efficiencies of MuSCAT. 


wavelength(nm) 

gW 

r^(%) 

Zs,2{%) 

wavelength(nm) 

gW 

r'(%) 

Zs,2{%) 

360 

0.00 

0.00 

0.00 

660 

0.00 

55.37 

0.00 

370 

0.00 

0.00 

0.00 

670 

0.00 

54.80 

0.00 

380 

0.00 

0.00 

0.00 

680 

0.00 

54.35 

0.00 

390 

0.00 

0.00 

0.00 

690 

0.00 

52.67 

0.00 

400 

2.02 

0.00 

0.00 

700 

0.00 

0.07 

0.00 

410 

42.21 

0.00 

0.00 

710 

0.00 

0.00 

0.00 

420 

44.63 

0.00 

0.00 

720 

0.00 

0.00 

0.00 

430 

46.75 

0.00 

0.00 

730 

0.00 

0.00 

0.00 

440 

48.39 

0.00 

0.00 

740 

0.00 

0.00 

0.00 

450 

50.36 

0.00 

0.00 

750 

0.00 

0.00 

0.00 

460 

50.99 

0.00 

0.00 

760 

0.00 

0.00 

0.00 

470 

52.42 

0.00 

0.00 

770 

0.00 

0.00 

0.00 

480 

52.35 

0.00 

0.00 

780 

0.00 

0.00 

0.00 

490 

53.23 

0.00 

0.00 

790 

0.00 

0.00 

0.00 

500 

53.56 

0.00 

0.00 

800 

0.00 

0.00 

0.00 

510 

54.21 

0.00 

0.00 

810 

0.00 

0.00 

0.00 

520 

54.45 

0.00 

0.00 

820 

0.00 

0.00 

0.08 

530 

54.87 

0.00 

0.00 

830 

0.00 

0.00 

35.05 

540 

54.58 

0.00 

0.00 

840 

0.00 

0.00 

36.12 

550 

15.81 

0.00 

0.00 

850 

0.00 

0.00 

34.50 

560 

0.00 

0.06 

0.00 

860 

0.00 

0.00 

32.74 

570 

0.00 

57.17 

0.00 

870 

0.00 

0.00 

30.78 

580 

0.00 

56.60 

0.00 

880 

0.00 

0.00 

29.08 

590 

0.00 

57.07 

0.00 

890 

0.00 

0.00 

27.13 

600 

0.00 

56.74 

0.00 

900 

0.00 

0.00 

25.21 

610 

0.00 

56.73 

0.00 

910 

0.00 

0.00 

22.99 

620 

0.00 

56.56 

0.00 

920 

0.00 

0.00 

14.40 

630 

0.00 

56.29 

0.00 

930 

0.00 

0.00 

0.02 

640 

0.00 

55.97 

0.00 

940 

0.00 

0.00 

0.00 

650 

0.00 

55.63 

0.00 

950 

0.00 

0.00 

0.00 


Table 3 Summary of linearity information of CCDs for each gain and each readout speed. 


Gain mode 

ADC speed [Hz] 

t/^-band 

linearity slope [%] 

r2-band 

inearity slope [%] 

2:s_2-band 

linearity slope [%] 

1 

2 X 10*^ 

0.154 ±0.062 

0.104 ±0.051 

0.204 ±0.026 


1 X 10^ 

-0.198 ±0.071 

-0.055 ±0.043 

-0.133 ±0.040 

2 

2 X 10® 

0.016 ±0.072 

0.018 ±0.045 

0.112 ±0.025 


1 X 10® 

-0.210 ±0.069 

0.005 ±0.043 

-0.061 ±0.039 

4 

2 X 10® 

0.169 ±0.044 

0.120 ±0.037 

0.115 ±0.030 


1 X 10® 

-0.059 ±0.077 

-0.085 ±0.043 

-0.037 ±0.034 
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Table 4 Summary of throughput (TP) of MuSCAT on the 188cm telescope at OAO 


band 

atmos. transmittance 

Ml 

M2 

Muscat 

expected total TP 

measured TP 

92 

~60%“ 

86%" 

79%" 

SlW 

21% 

20% 

r'2 

~65%“ 

85%" 

85%" 

56%'* 

27% 

28% 

4.2 

~80%“ 

~80%" 

~80%" 

29%'* 

15% 

13% 


Ml = main mirror, M2 = secondary mirror, MuSCAT = all optics * BEAR coating * QE, 


typical, measured, extrapolated from measured data up to 750 nm, expected 


Table 5 Summary of the test observation and analysis 


Target (Filter) 

Exp. time [s] 

^data ^ 

A^comp 

Ft/F, - 

^ap [pixel] 

rms ® [%] 

GJ436 (ff') 

30 

99 

2 

0.42 

24 

0.101 

GJ436 (r0 

30 

99 

3 

0.85 

26 

0.074 

GJ436 (z' 2 ) 

30 

100 

2 

3.6 

24 

0.076 

WASP-12 (5*) 

30 

79 

4 

0.38 

22 

0.16 

WASP-12(r^) 

30 

82 

3 

0.52 

24 

0.16 

WASP-12 (z', 2 ) 

60 

45 

2 

0.59 

22 

0.15 


^ The number of observed data points. 

^ The number of compaiison stars used for relative photometry. 

^ Unnormalized flux ratio of the target star and the ensemble of the comparison stars. 
^ Applied aperture radius. 

® The rms value of the residual light curve from a linear fit. 


Table 6 Error budget 


Target (Filter) 

'T’target [%] 

<7'comp [%] 

O-sky [%] 

^read [%] 

<^scin [%] 

'T'unknown [%] 

rms [%] 

WASP-12 {g' 2 ) 

0.058 

0.036 

0.092*^ 

0.005 

0.078 

0.076 

0.16 

WASP-12 {r' 2 ) 

0.047 

0.034 

0.077 

0.004 

0.081 

0.106*" 

0.16 

WASP-12 {zs,2) 

0.065 

0.048 

0.083*^ 

0.006 

0.057 

0.083*" 

0.15 

GJ436 {g' 2 ) 

0.043 

0.028 

0.037 

0.003 

0.050 

0.061*" 

0.101 

GJ436 {r' 2 ) 

0.023 

0.021 

0.019 

0.001 

0.050*" 

0.041 

0.074 

GJ436 (z,,2) 

0.023 

0.043 

0.023 

0.003 

0.050*" 

0.019 

0.076 


^ The bold text indicates the most dominant noise component for each light curve. 
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